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g '■ ABSTRACT 

o 

. We examine the cooling radiation from forming galaxies in hydrodynamic simula- 

tions of the LCDM model (cold dark matter with a cosmological constant), focusing 
O I on the Lya line luminosities of high-redshift systems. Primordial composition gas con- 

I denses within dark matter potential wells, forming objects with masses and sizes com- 

' parable to the luminous regions of observed galaxies. As expected, the energy radiated 

■ in this process is comparable to the gravitational binding energy of the baryons, and 

^ ' the total cooling luminosity of the galaxy population peaks at z ~ 2. However, in 

. contrast to the classical picture of gas cooling from the ~ 10^ K virial temperature 

of a typical dark matter halo, we find that most of the cooling radiation is emitted 
by gas with T < 20, 000 K. As a consequence, roughly 50% of this cooling radiation 
emerges in the Lya line. While a galaxy's cooling luminosity is usually smaller than the 
ionizing continuum luminosity of its young stars, the two are comparable in the most 
massive systems, and the cooling radiation is produced at larger radii, where the Lya 
photons are less likely to be extinguished by dust. We suggest, in particular, that cool- 
ing radiation could explain the two large (~ 100 kpc), luminous {Lj^ya ~ lO^^ergs"^) 
"blobs" of Lya emission found in Steidel et al.'s (1999) narrow band survey of a z = 3 
proto-cluster. Our simulations predict objects of the observed luminosity at about the 
right space density, and radiative transfer effects can account for the observed sizes and 
line widths. We discuss observable tests of this hypothesis for the nature of the Lya 
blobs, and we present predictions for the contribution of cooling radiation to the Lya 
luminosity function of galaxies as a function of redshift. 
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Subject headings: cosmology: theory — galaxies: formation — galaxies: ISM 

1. INTRODUCTION 

Radiative cooling of gas in dark matter potential wells is an essential element of the current 
theoretical understanding of galaxy formation. As the baryons condense into tightly bound clumps 
at the centers of these potential wells, they must radiate away the energy that they acquire through 
compression and shock heating. At first glance, it appears that this radiation must be negligible 
compared to the radiation produced by stars: the gravitational binding energy per unit mass is 
~ ~ lO^^c^ even for a (high) internal velocity v ~ 300 kms"-*^, while the nuclear energy per unit 
mass released in enriching gas to solar metallicity is ~ 10^'^ c^. However, for primordial composition 
gas cooling to T ~ 10*^ K, all of the emitted radiation is either hydrogen or helium line radiation 
(primarily hydrogen Lya) or continuum radiation above the hydrogen ionization threshold. Only 
the hottest stars emit photons at these energies. Furthermore, the cooling radiation is produced 
mainly at large radii, where it can escape the galaxy without being absorbed by dust. Finally, the 
kinetic energy of supernova explosions is deposited in dense, interstellar gas, where much of it can 
emerge in the form of hydrogen line emission and ionizing continuum radiation. 

Because of these factors, cooling radiation could make a significant contribution to the X-ray, 
UV, and hydrogen line emission from young galaxies. In this paper, we examine this contribution 
using hydrodynamic cosmological simulations. Although many such simulations incorporate radia- 
tive cooling, most studies have not followed the emitted radiation in detail. One notable exception 
is the work of Cen &; Ostriker (1992, 1996), who calculate the spectrum of the background produced 
by cooling radiation, but these papers do not examine the emission on a galaxy by galaxy basis. 
Closer to the focus of this paper is the work of Katz (1992), which examined the formation of a 
single, massive galaxy, finding that the amount of gravitational energy radiated in atomic lines was 
comparable to the amount of energy injected by supernovac. Here we examine the cooling radiation 
from the whole population of forming galaxies, though the individual objects in our simulations are 
not as well resolved as that of Katz (1992). 

The other main approach to theoretical modeling of galaxy formation utilizes semi-analytic 
methods, in the tradition of White &: Frenk (1991). These methods rely on simplifying assumptions 
about gravitational collapse and gas dynamics — in particular, they usually assume that infalling 
gas shock heats to the virial temperature of the dark matter halo, before cooling to join the central 
object. Our results below suggest that this assumption may break down in the messy assembly 
process that characterizes hierarchical galaxy formation. 

The observational study of star formation in high-redshift galaxies has accelerated in recent 
years, with progress in wavelength ranges from the UV/optical (Steidel et al. 1996; Williams et al. 
1996) to the IR (Elbaz 1999) to the sub-mm (Barger et al. 1999). The present paper is motivated 
mainly by recent progress in Lya emission line searches. Historically, this approach to finding high- 



-3- 



redshift galaxies has proved difficult (Pritchet 1994; Thompson, Djorgovski and Trauger 1995), 

in part because the Lya emission associated with star-forming regions may be damped by the 
combined effects of resonant scattering and dust absorption. However, Lya emission line searches 
have recently begun to bear fruit, thanks to improvements in sensitivity and efficiency (Hu, Cowie 
& McMahon 1998; Steidel et al. 1999; Rhoads et al. 2000). 

One of the most striking results of these searches is the discovery by Steidel et al. (1999, 
hereafter S99) of two large regions of diffuse Lya emission, in their narrow band observations of a 
galaxy proto-cluster at z = 3.1. These "blobs" do not resemble the typical galaxies detected in the 
proto-cluster. They have characteristic angular sizes of 15" and fluxes of 1.3 x 10~^^ ergcm"^ s~^. 
For the cosmological model that we adopt in this paper (fJ^ = 0.4, JIa = 0.6), these properties 
imply proper linear sizes of I ~ 75/i~^kpc and luminosities of L^ya ~ 4.4 x 10^^ ergs~^.^ 
Some other workers may have observed similar objects. Francis (1999) found diffuse Lya blobs in 
a cluster at z = 2.4. Keel et al. (1999) found three diffuse regions of Lya emission with sizes of 
~40/i~^kpc and luminosities of ~ 10^^ ergs""^ (in our adopted cosmology), also at z = 2.4. 
These are somewhat smaller and less luminous than the blobs found by S99, and some are located 
around known AGN, so they may or may not be related phenomena. 

Several explanations for these Lya blobs have been proposed. They could be caused by pho- 
toionization from a hidden central AGN (S99). High-redshift radio galaxies often show emission 
line regions with sizes and luminosities comparable to the S99 blobs, and the absence of radio 
lobes could simply indicate a radio-quiet AGN; the anisotropic emission of the standard unified 
AGN model would then provide an explanation for the absence of nuclear optical emission. An- 
other possibility, also discussed in S99, is that the gas in the blobs is photoionized by ultraviolet 
(UV) continuum radiation from young stars. The required star formation rate is ~30 f~l Mq yr~^, 
where /e^c is the fraction of hydrogen ionizing radiation that escapes from the galaxy. Low-redshift 
observations imply upper limits of fesc ^ 10% (Leitherer et al. 1995; Hurwitz, Jelinsky &: Dixon 
1997), so unless these presumably gas-rich galaxies have much larger /esc, the required star forma- 
tion rates are substantial. However, current limits are quite weak, allowing rates < lOOOM0yr~^ 
(Steidel et al. 1999). Taniguchi &: Shioya (2000) have suggested that the Lya emission emerges 
from a galactic superwind that is powered by high star formation rates and extends to scales that 
are even larger than those of winds observed at low redshift (Heckman, Lehnert, & Armus 1993). 
While the luminosities and sizes of these blobs can be explained by wind models, the mechanism 
for channeling the wind energy into Lya emission is somewhat unclear. 

Drawing on the results of our more general investigation, we propose an alternative explanation 
for the Lya blobs: they represent gas that is radiating away its gravitational potential energy as 
it settles into massive galax;ies. We describe our numerical simulation techniques, with particular 
attention to our treatment of radiative cooling, in §2. In §3 we present our results for the cooling 
radiation and Lya emission from forming galaxies. In §4 we discuss whether the cooling radiation 



^With Q,m = 1, f^A = 0, these numbers become I ^ 55h ^ kpc and LLycx « 2.4 x 10''^ h ^ ergs 
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model can account for the observed properties of the Lya blobs. Section 5 discusses some of 

the numerical uncertainties in our calculations, compares our model to a similar model proposed 
independently by Haiman, Spaans, & Quatacrt (2000) using semi-analytic methods, summarizes 
our results and discusses their implications and predictions. 

2. SIMULATIONS 

We perform our simulations using the parallel version of the cosmological N-body/hydrodynamic 
code TreeSPH (Hernquist & Katz 1989; Katz, Weinberg, & Hernquist 1996a; Dave, Dubinski, & 
Hernquist 1997), a code that unites smoothed particle hydrodynamics (SPH) (Lucy 1977; Gingold 
& Monaghan 1977) with the hierarchical tree method for computing gravitational forces (Barnes 

& Hut 1986; Hernquist 1987). Dark matter, stars, and gas are all represented by particles; colli- 
sionless material is influenced only by gravity, while gas is subject to gravitational forces, pressure 
gradients, and shocks. The gas can also cool both radiatively, assuming primordial abundances, 
and through Compton cooling. 

TreeSPH is fully adaptive in both space and time. In SPH, gas properties are computed 
by averaging or "smoothing" over a fixed number of neighboring particles, 32 in the calculations 
here. Hence the smoothing lengths in TreeSPH decrease in collapsing regions, in proportion to 
the interparticle separation, and in underdense regions the smoothing lengths are larger. TreeSPH 
allows particles to have individual time steps according to their physical state, so that the pace 
of the overall computation is not driven by the small fraction of particles requiring the smallest 
time steps. There is a maximum allowed timestep, called the system timestep, and all particles are 
integrated with this step or one a power of two smaller. The timestep criteria are detailed further 
in Katz, Weinberg, &; Hernquist (1996a) and Quinn et al. (2000); we set the tolerance parameter rj 
to 0.4. 

Since this paper concerns cooling radiation, it is important to understand how we evolve the 
thermal energy. Usually one has two choices: to integrate the thermal energy equation explicitly 
or implicitly. Integrating explicitly would require a timestep about three times smaller than the 
cooling time, which is much smaller than what would be required to integrate the dynamical equa- 
tions, making the calculation prohibitively expensive computationally. To integrate the equations 
implicitly involves inverting an Ngas x Ngas matrix several times per timestep. With Ngas > 3 x 10^ 
in our simulations, this would also be prohibitively expensive. Instead we take an intermediate 
approach, solving the thermal energy equation scmi-implicitly as described in Hernquist and Katz 
(1989). Briefly, we integrate the changes in the thermal energy caused by shocks and pressure 
forces explicitly, while we integrate those caused by radiative cooling implicitly. Since the radiative 
cooling depends only on a particle's own temperature and density, independent of other particles, 
no matrix inversion is required, and the time scale of the non-radiative processes is comparable to 
the other time scales that govern the gas timestep. 



-5- 



The semi-implicit approach guarantees that we integrate the thermal energy equation in a 

stable way, but it does not guarantee accuracy. Accuracy is maintained by keeping the size of the 
timestep reasonably small. We accomplish this in two ways. First, we integrate the thermal energy 
equation for all gas particles, independent of their dynamical timestep, using a timestep that is one 
half the size of the smallest dynamical timestep of any particle. Second, we damp the cooling rate 
so that no gas particle loses more than a given fraction of its thermal energy in one timestep (see 
Katz & Gunn 1991); i.e., wc slightly change the physics to make the numerical integration more 
robust. This could have the effect of making some regions temporarily hotter than they would have 
been if these numerical compromises were not made. In practice, it just makes some regions take 
two or three timesteps to cool instead of one, since these regions have such short cooling times. 
This could change the temperature at which the energy is radiated, but it should not change the 
total radiated energy by a large amount. 

The simulations calculate radiative cooling processes assuming a primordial composition gas 
and ionization equilibrium. These processes include collisional excitation, radiative and dielectronic 
recombination, collisional ionization, and bremsstrahlung. Molecular processes and metal lines are 
omitted, which means that radiative cooling alone cannot reduce the temperature below ~ 10^ K. 
The gas cooling rates we use for these processes are listed in Katz, Weinberg, & Hernquist (1996a, 
see their figure 1). However, the rate at which the gas radiates energy is somewhat different from 
the rate at which it loses thermal energy: the emitted recombination radiation includes both the 
actual gas cooling and the atomic ionization potential, while collisional ionization removes thermal 
energy but yields no radiation at all. A possible deficiency of the simulations is the omission of 
He I line cooling, but we have verified that it is unimportant compared to the other processes in 
the simulations presented here. 

For the purposes of this paper, we are particularly interested in the Lya emission. We assume 
that a fraction 0.68 of recombinations to H I produce a Lya photon, appropriate for optically thick 
gas at 10^ K; accounting for the relative photon energies this channels 0.49 of the recombination 
energy through Lya. In practice this is usually a small contribution compared to the collisional 
excitation. Excitation of H I by collisions with electrons can result in either a Lya photon or in 
2-photon decay from the mctastable 2 state, as well as photons from higher series in some cases. 
The ratio between Lya and 2-photon decay is weakly dependent on temperature. By averaging over 
one simulation, we find that a fraction 0.59 of the H I collisional excitation is channeled through 
Lya with almost all of the remaining energy in the 2-photon continuum. 

TreeSPH can also include a metagalactic ionizing UV background field. TreeSPH does not 
include radiative transfer, so even dense regions, which in reality should be self-shielded, are exposed 
to the full background field. This results in a large energy exchange and unphysical gas radiation 
from these regions. Hence, we have restricted our analysis to simulations performed without such 
a background, which should yield reliable results for our purposes. 

In these simulations we include star formation and supernova feedback using the algorithm 
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described in Katz, Weinberg, k, Hernquist (1996a). The algorithm forms stars in dense cool regions 

at a rate essentially controlled by the dynamical supply of gas. For each star formation event, 
supernova energy is added to the surrounding particles as thermal energy on a time scale of 2 x 10^ 
years. 

In the version of our code used for these simulations, the thermal energy input from supernovae 
is added at the beginning of each system timestep, not continuously. The system timesteps are 
always larger than the cooling time in these dense star forming groups, so the excess thermal energy 
has radiated away by the end of the step when we output the system state. This actually turns out 
to be quite helpful, since the remaining radiation is caused by the radiating away of gravitational 
energy and thus can be measured in isolation. We can easily recover the supernova thermal input 
because it is proportional to the star formation rate. Using a Miller-Scalo initial mass function 
(IMF) with cutoffs at 0.1 and 100 M©, the heat input from supernovae is 2.7 x 10^^ ergs~^MQ^yr 
when smoothed over the neighboring particles (the exact value would be 2.5 x 10^^ ergs~^MQ^yr 
if there were no smoothing). The drawback to this method is that we cannot easily compute the 
distribution of the supernova energy within galaxies or the processes by which it is radiated away. 
Since the supernova energy is deposited in the dense ISM over a length scale determined by our 
resolution, its space and temperature distributions are somewhat suspect in any case, and we will 
not attempt to reconstruct them. 

In addition to cooling radiation, we compute the photoionizing radiation from hot stars from 
the simulation outputs. Using the code STARBURST99 (Leitherer et al. 1999), we find a Lyman 
continuum intensity of 2.6 x 10^^ ergs^^Mg^ yr for a Miller-Scalo IMF, assuming a mean Lyman 
continuum photon energy of 1.4 Ryd. For gas at 10^ K, there are 0.68 Lya photons emitted per pho- 
toionization (Chariot k. Fall 1993), giving an associated Lya emissivity of 9.4 x 10^^ ergs~^MQ^ yr. 
Because the stars that produce photoionizing radiation are short-lived, we can take the stellar 
emissivity to be proportional to the instantaneous star formation rate, which we compute from the 
gas distribution. We do not include this UV emission or stellar winds as a source of feedback in 
the simulations. It is likely to have even less dynamical impact than supernova feedback, since 
any energy captured by the gas is deposited in the densest, star forming regions with moderate 
temperatures, where it can be quickly radiated away. 

To summarize, the three sources of radiative energy included in these simulations that could 
result in Lya radiation are gravitational cooling, supernova feedback, and photoionizing input from 

hot stars. These sources have distinct physical origins, our code follows them separately, and we 
will keep them separate in the discussions and figures. When we refer to the "total cooling", we 
mean the gravitational cooling summed over all radiative processes, not that we are including the 
supernova or photoionizing input. 

We calculate the cooling radiation at discrete output times with the analysis program TIPSY,^ 



^http://www-hpcc. astro.washington.edu/TSEGA/tools/tipsy/tipsy.html 
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using the same cooling algorithm as TreeSPH. To identify discrete objects in the simulations, we use 
the program SKID.^ Because SKID defines gravitationally bound groups of particles, we will refer 
to these objects as "groups," though each such group actually represents a single galaxy. SKID 
slides particles meeting density and temperature criteria along density gradients until they reach a 
local maximum. We use a density cutoff corresponding to the edge of a virialized isothermal halo, 
or a gas density of pg = {ilb/^m){Pvir/^)j where the mean density of a spherical virialized halo 
is 178 at early times and is given in general by Kitayama &: Suto (1996). As in Weinberg et al. 
(1999), we also restrict ourselves to groups that contain more than 64 baryonic particles, roughly 
the mass resolution limit of the simulations. Once we identify the groups, we can easily add up the 
cooling radiation for all particles in the group. Cooling radiation can also be emitted outside these 
groups, but it proves to be a minor contribution except at very early times; most of this additional 
emission comes from groups excluded by our mass resolution criterion rather than from low density 
gas. 

All three of the simulations we discuss in this paper assume a A-dominated cold dark mat- 
ter cosmological model with = 0.4, f^A = 0.6, h = iJo/(100kms~^ Mpc~^) = 0.65, and a 
primeval spectral index n = 0.93. With the tensor mode contribution, normalizing to COBE using 
CMBFAST (Scljak & Zaldarriaga 1996; Zaldarriaga, Seljak, & Bcrtschingcr 1998), implies a nor- 
malization (Tg = 0.8, which provides a good match to cluster abundances (White, Efstathiou, & 
Prenk 1993). We use the Hu & Sugiyama (1996, equation D28) formulation of the transfer func- 
tion. We adopt a baryonic density = 0.02 h"^ consistent with the deuterium abundance in high 
redshift Lyman limit systems (Buries & Tytler 1997, 1998) and with the opacity of the Lya forest 
(Ranch et al. 1997). All of our simulations model a triply periodic cubical volume. 

The main and largest simulation we discuss is the L144 simulation (Dave et al. 2000), with 144^ 

gas and dark matter particles, a box length of 50 comoving Mpc on a side, and a gravitational 
softening length €grav = 7 h^^ comoving kpc (equivalent Plummcr softening). The nominal gas 
mass resolution is 5.4 x 10^*^ Mq, corresponding to 64 gas particles. To investigate the effects of 
our finite resolution we perform two additional simulations. The Lll/64 simulation has a higher 
spatial and mass resolution, so it must be run in a smaller volume, 11.1 comoving Mpc on 
a side. It has 64^ gas and dark matter particles, egrav = 3.5 comoving kpc, and a gas mass 
resolution of 6.8 x 10^ Mq. The Lll/32 simulation is identical to the Lll/64 simulation, with the 
same initial phases, except that it uses 32^ particles of each type and has the same resolution 
as the L144 simulation. In all the simulations the nominal spatial resolution in physical units is 
~ 2egrav/(l + z). 



^http://www-hpcc. astro.washington.edu/TSEGA/tools/skid. html 
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3. COOLING RADIATION AND THE LYa LUMINOSITY FUNCTION 

Figure 1 displays the temperature distribution of the total gravitational cooling radiation (all 
radiative processes), in the L144 simulation at z = 3 and z = 0. In the conventional theoretical 
sketch of galaxy formation (e.g., White &; Rees 1978; White &; Prenk 1991), gas that falls into dark 
matter halos is shock heated to the virial tcmpcratTirc T^,j,. = {limp/2k)vl « 106K(vc/165kms-^)2, 
then cools and settles into the central galaxy. The dotted line in Figure 1 shows the temperature 
distribution of cooling radiation expected for gas cooling from an initial temperature of 10^ K. The 
numerical simulation results paint a very different picture. A large fraction of the cooling radiation, 
75% at z = 3 and 30% at z = 0, comes from gas with 10^ K < T < 2 x 10^ K. Most of the remaining 
radiation comes from much hotter gas, with 10^ K < T < 10^ K. 

The lack of cooling radiation from gas with 2 x 10^ K < T < 10^ K cannot be a result of gas 
"moving quickly" across this temperature range, since even if it did, it would still have to radiate 
away its thermal energy in order to reach 2 x 10^ K. Instead, Figure 1 implies that most of the gas 
that cools into galaxies is never heated to the virial temperature of a galaxy-mass dark halo. This 
conclusion accords with that of Kay et al. (2000), who find, based on similar sorts of simulations, 
that only ~ 10% of SPH particles that end up in galaxies were ever heated above 10^ K. We find 
qualitatively similar results in our own simulations, but we have not examined particle temperature 
trajectories with the high (single timestep) time resolution used by Kay et al. (2000). 

This physical result has major implications for the spectrum of cooling radiation from forming 
galaxies, since much of the energy from neutral gas at T ~ 10^ K emerges in the Lya line as a 
result of collisional excitation (see Katz, Weinberg, & Hernquist 1996a, figure 1). The hot gas 
cooling, on the other hand, is dominated by bremsstrahlung. At z = 3, Lya emission accounts 
for 43% of the cooling radiation in the simulation, with 29% emerging in HI 2-photon emission, 
19% in bremsstrahlung, and only 8% in other radiative processes. If we assumed a metal-enriched 
intergalactic medium instead of primordial composition, then more of the hot gas might be able to 
cool, boosting the fraction of radiation from 10^ K < T < 10^ K, but the large emissivity of lower 
temperature gas would remain. 

Our focus in this paper is the cooling radiation, and especially the Lya emission, associated 
with individual galaxies. In Figure 2, the dot-dashed line shows the cumulative luminosity function 

of the total gravitational cooling at z = 3. By "total" we mean that this includes all of the cooling 
radiation processes, not just Lya radiation. The reradiated supernova energy is plotted as the 
dotted line. The gravitational energy available to a group increases more quickly with mass than 
the star formation rate. The dominant source of the total cooling radiation thus changes over from 
reradiated supernova energy to gravitational energy as the mass and luminosity increase. The two 
processes combined give the luminosity function shown by the solid line. Finally, the dashed line 
shows the stellar UV emissivity, calculated from the star formation rate in the galaxy. For all 
luminosities, there is more energy in the UV radiation produced by the massive stars than in the 
cooling radiation. However, the UV radiation from young stars is likely to be heavily absorbed by 
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dust, and reradiated in the infrared. The gravitational coohng radiation emerges from lower density 

gas and is more likely to escape. The supernova energy may also stand a much better chance of 
escaping, as supernovae can destroy the dense clouds responsible for the heaviest absorption and 
deposit their kinetic energy in a more diffuse medium. 

We plot the cumulative luminosity function of the Lya cooling radiation alone in Figure 3, 
at several redshifts. Since the fraction of supernova energy converted into Lya is uncertain, we 
include only gravitational sources of cooling in this plot. Figure 2 shows that gravitational cooling 
dominates in the most luminous objects. The number of highly luminous objects reaches a peak 
at z = 2, and declines thereafter. At ^; = 3 there are ~4 x 10~^ Mpc~^ objects in the simulation 
with Lya luminosities greater than 3 x 10^^ /i~^ergs~^, comparable to the 'blobs" of S99. 

In the left hand panel of Figure 4, we show a map of the Lya emission contributed by grav- 
itational energy in one of our simulated groups, and in the right hand panel we plot the stel- 
lar photoionization using a conversion factor of star formation to recombination-induced Lya of 
6.6 X 10^^ ergs~^MQ^ yr. The emission is shown by the gray scale images, with the intensity scale 
marked in terms of logj^ol-^ijya/l^'^gs"^ cm~^ asec~^)]. The emission from gravitational cooling is 
more spatially extended than the stellar emission. Emission of cooling radiation from reradiated 
supernova feedback would look like a scaled version of the right hand panel, slightly smoothed by 
the feedback algorithm. We also plot contours of the neutral hydrogen column density, calculated 
assuming an ionizing background of 3 x 10^^ ergs~^ cm~^ Hz~^ sr~^ and correcting for self-shielding 
(Katz et al. 1996b). Self-shielding accounts for the rather sharp edges seen in the neutral gas (cf. 
Maloney 1993). These emission plots ignore radiative transfer effects and dust extinction, both of 
which could greatly alter the observed appearance of these systems. We discuss these issues in §4. 

The spatial extent of the emission for the groups in the L144 simulation is shown in Figure 5. 
The measure we use here is the root mean square distance from the center of the group, weighted 
by the Lya emissivity. The dotted line shows the gravitational softening length; recall that the 
nominal spatial resolution is about twice this length. We plot the Lya extent against the rnis extent 
of the gas mass defined in a similar manner. We also show the size weighted by star formation rate, 
which would represent the size of the stellar UV emissivity region. In all cases, the star formation 
in our simulations occurs in a small, partially resolved region at the center of the group. The Lya 
emission usually emerges over a region comparable to the size of the gas as a whole and is well 
resolved numerically. However, the most luminous groups in the simulation tend to have more 
concentrated emission; although there is weak emission at large radii, the typical rms sizes are less 
than 10 kpc. This size is much smaller than that of the blobs observed by S99, a point we will 
return to in §4. 

We plot the gravitational cooling in Lya as a function of galactic mass in Figure 6. We 
define the galactic mass to be the stellar mass plus the mass of the gas that is at least 1000 times 
overdense and has temperature T < 3 x 10^ K. To increase our dynamic range, we show both 
the L144 simulation and the higher resolution (but smaller volume) LI 1/64 simulation. There 
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is a strong correlation between the emitted Lja cooling radiation and the galactic mass. In the 
L144 simulation, the galaxies with high Lya cooling luminosities are all high mass objects, many 
corresponding to galaxies or above. 

In Figure 7 we plot the total gravitational cooling radiation as a function of the star formation 
rate at z = 3. Supernova feedback would add an amount of cooling radiation shown by the dashed 
line. The Lyman continuum luminosity associated with the stellar UV emission is marked along 
the top axis. Star formation rates are underestimated in marginally resolved systems, a numerical 
artifact that causes the spread in these scatter plots towards low star formation rates at relatively 
low cooling luminosity. Figure 7 demonstrates on an object-by-object basis the features seen in the 
luminosity functions of Figure 2. Supernova cooling dominates gravitational cooling in low mass 
objects, but gravitational cooling takes over at high masses. The objects with the highest cooling 
luminosity also have high star formation rates (SFR > lOOMQyr^^). The stellar UV luminosity 
always exceeds the gravitational cooling radiation (all points representing well resolved objects lie 
below the dotted line), but they are of similar magnitude in the most luminous objects, so the 
cooling radiation would dominate if the stellar UV is heavily extinguished by dust. 

To examine the energetics of the cooling radiation, we would like to define the gravitational 
potential energy available to the baryons in the groups, but this definition is quite ambiguous. 

Rather than being symmetric, isolated entities, the groups are embedded in a complex, filamentary, 
and clumpy structure. We are interested in the energy available to the baryons, but they interact 
with the dark matter in a complex manner and energy is transferred between them. Nevertheless, 
it is interesting to see whether the potential energy, with the matter outside the group ignored, is 
a good predictor of the cooling radiation. 

In all of the previous figures, the groups defined by SKID consist of baryonic particles only. 
To study the gravitational energy, we must redefine the groups to include dark matter particles as 
well, using a density cutoff of ptot = Pvir/^- If we defined Ufot to be the total potential energy of the 
group, and if the baryons and dark matter were distributed identically with a baryon fraction /;,, 
and if we assigned half of the baryon-dark matter interaction energy to the baryon potential energy 
Ub, then we would have Uh = fbUtot- Another possibility would be to assign the entire interaction 
energy to the baryons. For example, if the cooling and collapse of baryons to the center of the 
group takes place when the dark matter halo is already assembled, then this latter definition is 
closer to the amount of energy that the gas radiates. As a practical matter, M^V^, where Vc is the 
circular velocity and Mj, is the baryonic mass, usually lies between these two definitions. Hence, 
we take MhV^ as our estimate of the potential energy of the baryons. One might object that, by 
the virial theorem, the energy available is only \Uh\/2. However, the virial theorem does not apply 
to the baryons alone, as they are confined by the potential of the dark matter.^ As a rule the sum 
of the kinetic and thermal energies falls short of even the smaller definition of |?7b| by a factor of 2 



*In addition, the group as a whole is not an isolated system, and in our simulations the gravitational force is 
softened and thus not a power law; both of these facts also violate the conditions for the virial theorem to hold. 
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to 4. 

A simple estimate of the cooling radiation is then Mi,V^/t, the potential energy divided by the 
cosmic time. At z = 3 this prescription somewhat underestimates the cooling from each group, as 
shown in Figure 8, though the discrepancy is not large considering the ambiguity in our definition 
of \Ub\. At z = 0, this simple formula overestimates the cooling, indicating that the cooling occurs 
preferentially at early times. There could be several reasons for this trend. The cooling could be 
taking place less effectively since the virial temperatures are rising; the rate of accretion of matter 
into the group could have slowed; or the baryons could become more effective at transferring their 
energy to the dark halos. We are unable to determine the answer from the discrete outputs of the 
current simulations. Tracking the evolution of distinct groups as a function of time, as was done 
by Katz (1992) for a single galaxy, would give more insight into the energetics, but we leave this 
investigation for future work. 

4. MODELING THE LYa BLOBS 

So far, we have established that galaxy-like groups in our simulations generically show large 
amounts of cooling in Lya. Can this cooling explain the "blobs" observed by S99? 

The Lya luminosities of some of our groups are as large as those of the S99 blobs. The 
number density of groups in our simulations with a Lya luminosity of 3 x 10^^/i~^ ergs~^ is about 
4 X 10~^ Mpc~^. At present we can make only crude estimates of the number density of the 
observed Lya blobs. The volume examined by S99 has a size of 8'.7x 8'.9 and a depth of = 0.066, 
or a comoving volume of 4200/t~^ Mpc^ in our adopted cosmology. As long as the blobs are less 
abundant in other regions of space, they have a comoving density of < 5 x lO^^/i^ Mpc"'^. A better 
estimate can be obtained from S99's observation that this region is overabundant in both Lyman 
break galaxies and in Lya-emitting galaxies by a factor of ~ 6. If the blobs are biased in the 
same way, they have a density of ~ 8 x 10~^^^ Mpc~^, twice the density we find. If these blobs 
are associated with very massive objects, they are probably more highly biased than Lyman break 
galaxies; hence their actual number density is probably lower than this estimate. With only two 
objects, which are probably correlated, the statistical uncertainties are large. So while the current 
constraints are quite weak, the number density is consistent with that found in our simulations. 

The most luminous objects in our simulations typically have high mass and are strongly clus- 
tered. For example, the 100 most massive objects in the L144 simulation have a comoving corre- 
lation length of 5.7 Mpc. This strong clustering is consistent with the discovery of the blobs in a 
proto-cluster, though of course the abundance of blobs in blind fields is not well known. 

As shown in Figures 4 and 5, the length scale over which most of the Lya cooling is emitted 
in our simulations is usually much smaller than the observed sizes of the blobs in S99, though it 
is also larger than the typical size of the star forming regions. Hence, we must appeal to resonant 
scattering to transport the Lya photons to large radii. 
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The radiative transfer of Lya in a static slab is a classic problem (Adams 1972, 1975; Neufeld 
1990). The behavior of Lya photons depends upon the optical depth of the slab. For intermediate 
optical depths, the Lya photons are spatially trapped, and they escape by scattering into the tails of 
the Doppler distribution, where the slab is optically thin. For very large optical depths, the photons 
escape by scattering into the damping wings, where they then perform a random walk in both space 
and frequency. The transition between these two regimes occurs roughly where the damping wings 
become optically thick. For gas with a temperature T4 = Tj (10^ K), the total central optical depth 
is TO = {Nm/l.^ X 10^3 cm-2)T4"^''^ and the damping constant is a = 4.7 x IQ-'^T^''''^ . The 
damping wings, which have the line profile 0(a;) = a/{nx'^), become optically thick at tq ~ 4 x 10^, 
or Nhi « 7 X cm'^. 5 

Despite the extensive work on the slab problem, we have not found calculations of the typical 
line-center optical depth at last scattering (riast) in the literature. We can estimate Tiast in two 
regimes. In the case of intermediate optical depth, the spatial diffusion is negligible and Tiast ^ tq/2, 
the optical depth to the slab center. The photons escape with a double-peaked profile with a 
typical frequency shift x = Aiy/AuDop = V^nro. For large optical depths, the random excursion 
that leads to escape takes place at a frequency ~ (aro/2-^/7r)^/^ (Adams 1972), for which the 
optical depth through the slab is r* = To0(a;*) ~ (aro/2-^/7r)^/'^. The last scattering of the photons 
in that excursion occurs at optical depth ~ 1 at that frequency, or a line-center optical depth 
of Tiast ~ (t>ix*)^^ ^ (ro/2)(aro/2A/7r)^^/^. Without performing Monte Carlo calculations, it is 
difficult to say exactly where the transition between these regimes occurs; we will take it to be 
To ~ 10^. For very large columns, the Lya photons arc extinguished by dust. This occurs for 
^Hl ^ 4 X 10^*^ ^^~'^'^4^'^^dusi^ where S,dust is the dust-to-gas ratio relative to the Galactic value 
(Ncufcld 1990). Even though galaxies at high rcdshift might have low metallicity, star-forming 
regions should still be fairly rich in dust. We somewhat arbitrarily take £,dust = 0.1 as a typical 
value. 

Most of the neutral gas in our baryonic groups is at about 10^ K. Hence, we might expect 
that Lya photons emitted at iV > 2 x 10^^ cm~^ are absorbed by dust. Photons emitted at 
2 X 10^^ cm~^ < Nhi ^ 2 x 10^-^ cm~^ are scattered in the line wings and finally escape in the 
range 7 x 10^^ cm~^ < Nhi ^ 7 x IQ-*^^ cm~^. Finally, photons emitted at Nhi ^ 2 x 10^^ cm~^ are 
scattered mostly in place and eventually escape not far from their region of formation. 

Although we are using one-dimensional models for these estimates, they should be at least 
crudely applicable to our three-dimensional groups. Applying this picture to Figure 4, we see that 
Lya photons originating from stellar photoionization are likely to be absorbed. This is consistent 
with the strong absorption seen in star-forming galaxies at high redshifts (Small, Ivison & Blain 
1997). However, many of the photons from gravitational cooling may be scattered outwards to the 
lO^^cm"^ neutral hydrogen contour or about ~30/t~^kpc, with velocity widths up to 300kms~^. 



""Note that this column density is much less than what is conventionally referred to as a "damped" line, Nhi ~ 
102°cm-2. 
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In general we find that much of the gravitational cooling is emitted outside the region of intense 
star formation. Usually ~ 20% is generated outside our nominal dust cutoff, but this fraction can 
be smaller or larger depending on the value wc choose for the cutoff. A crucial ingredient in this 
argument is the presence of reservoirs of neutral gas at large radii around the largest groups in our 
simulations. The velocities obtained are not quite as high as in Figure 8 of S99, but the greatest 
velocity spread there comes from isolated knots and the velocity of the truly diffuse emission is 
mostly unconstrained. 

Clumping of the gas could make the escape of the photons easier and shrink the apparent size 
of the Lya-emitting regions (Neufeld 1991). The presence of large bulk motions can greatly affect 
the line transfer, probably explaining the moderate fraction of cases where photons escape from 
star-forming regions (Kunth et al. 1998; Ahn &: Lee 1998). However, calculations using the velocity 
fields of our groups, which are not well resolved in any case, are beyond the scope of this paper. 
We can conclude at this point only that our model may be consistent with the observed sizes and 
velocity widths of the Lja blobs. 

The cooling gas model for the blobs makes several testable predictions. Since the Lja emis- 
sion is expected to be scattered to larger radii, emission in other lines should be more centrally 
concentrated. The neutral hydrogen column should be ~ 10^^ cm~^ out to the radius of the Lja 
blobs. Since the Lya is caused mostly by collisional excitation rather than photoionization, the Ha 
flux should be quite small relative to Lya (< 2%). These predictions are at odds with the results 
of Francis et al. (1996) and Francis, Woodgate & Danks (1997), who report a detection of Ha and 
CIV in their three blobs and He H Balmer-a in two. In their brightest blob, labeled Bl, Ha has 
a similar strength and distribution to Lya. This suggests that at least these blobs are not due to 
gravitational cooling; this is concordant with the apparent double-lobed morphology of Bl and red 
stellar colors of the central galaxies, which suggest a AGN origin for these blobs. Ha Tinfortunatcly 
falls in the K-band only in the range 2.0 < z < 2.6, so it is important to search for other possible 
lines in the S99 blobs. 

Another implication of our model is that large luminosities in diffuse Lya emission are associ- 
ated with massive objects with large star formation rates. In fact, if Figure 7 is taken at face value, 
the star formation rates implied for the S99 blobs are ~lOOM0yr~^. However, the tight relation 
between cooling luminosity and star formation rate is due to a similarly tight relation between 
baryonic mass and star formation rate. These quantities may not be as correlated in real galaxies. 
The more fundamental prediction is that there is a massive galaxy at the heart of each Lya blob. 

As mentioned before, to power the Lya nebulosity with Lyman continuum radiation from 
young stars, one needs ~ 300 (/esc/0. 1)"^ yr'^ where / 

esc is poorly known. In the supernova 
wind model of Taniguchi & Shioya (2000) , the required star formation rates arc ~ 200 Mq yr~^ . This 
rough agreement between the required star formation rates of three different models is somewhat 
frustrating, with only the model of AGN photoionization allowing a different rate. 
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5. DISCUSSION 

The most basic, and perhaps the most surprising result of this work, is that most of the 
gravitational cooling radiation comes from gas at T < 2 x 10^ K, far below the typical virial 
temperatures of galaxies (see Figure 1). In these simulations (and those of Kay et al. 2000), most 
of the gas that settles into galaxies is never heated to the virial temperature of a galaxy-mass dark 
matter halo. In §5.1 we discuss some of the numerical issues surrounding this result. Astronomical 
implications are discussed in §5.2. 

5.1. Numerical Uncertainties 

While numerical simulations arc a useful guide to physical intuition, they are far from a perfect 
model of reality. With regard to our principal result — the temperature distribution of the cooling — 
it is easy to think of some limitations of the current code that might affect cooling estimates. For 
example, the omission of molecules and metals might have some effect. The main effect of molecular 
cooling would simply be to allow cooling below 10*^ K in the densest regions, so it would not affect 
our Lya predictions. Metals might channel some of the Lya cooling into other lines, and it would 
allow more high temperature gas to cool. However, the mctallicity must be above a few percent of 
solar before it significantly affects the cooling curve, and such a high metallicity seems somewhat 
unlikely for gas that is falling into galaxies for the first time at high redshift. The key feature that 
leads to high Lya luminosity is that most gas that enters high-redshift galaxies is never heated 
to high (T > 10^ K) temperatures at all. In this regard, a more worrisome concern is inadequate 
resolution of shocks, which could allow gas temperatures to rise more slowly than they should. 

Concern about numerical resolution in general can be addressed by examining the simulations 
Lll/32 and Lll/64. As we mentioned earlier, Lll/32 has the same resolution as our main simu- 
lation L144, but it covers a smaller volume. Lll/64 resimulates the same volume as Lll/32 with 
the same initial density field, but with eight times the mass resolution and twice the spatial resolu- 
tion. These two resolutions are compared in Figure 9, where we plot the luminosity of the cooling 
radiation emitted in Lya for both simulations. Since the initial density fields are identical, we can 
compare the Lya cooling luminosity galaxy by galaxy. We find that although the individual Lya 
cooling luminosities can vary, the ratio of Lya cooling luminosity in the two simulations scatters 
around one, giving us confidence in our ability to predict the Lya cooling luminosity in the L144 
simulation. The amount of scatter is not surprising in view of the highly stochastic nature we find 
for the cooling luminosity in individual simulated galaxies. 

Though it is encouraging that the cooling luminosity function is consistent between our two 
different resolutions, the cooling could be altered by resolution effects if the important scales are 
below the resolution of even the Lll/64 simulation. This caveat may undermine our prediction that 
Lya should dominate Ha and X-ray emission, as the presence of stronger shocks would increase 
collisional ionization and recombination radiation. We will eventually be able to test for resolution 
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effects over a wider dynamic range, but this will require other, more computationally demanding 
simulations. 

Our method for time integration of the thermal energy also deserves comment. As mentioned 
above, we damp the cooling when the cooling times are shorter than the minimum timestep. It 
would be surprising if this substantially affects the total radiated energy. However, if there is an 
induced error, it is in the direction of causing us to slightly overestimate the temperature of the 
gas and hence slightly overestimate the importance of X-ray emission relative to Lya emission. 

5.2. Observational Implications 

Laying aside our numerical concerns, which can only be fully addressed by future work, let 
us summarize our principal astronomical results. We find that large amounts of gravitationally 
induced Lya radiation should be produced around massive galaxies in the early universe. The 
gravitational cooling radiation is smaller than the stellar UV output, but its relative significance 
increases with the galaxy mass, and its physical extent is quite different, making it less likely to be 
extinguished by dust. The Lya luminosity function extends up to ~10^^ergs~^, and it reaches a 
peak at z ^ 2. 

We suggest that this gravitational cooling radiation could explain the blobs found by S99. The 
luminosities, number densities, and clustering are consistent with current constraints. The cooling 
radiation is dominated by coUisional excitation, so that the expected ratio of Ha to Lya is small. 
The size of the Lya emission region is unlikely to be as large as the observed blobs, but we find 
that sizes approaching those observed can be produced by resonant scattering of the Lya radiation. 

If our model of the blobs is correct, it follows that the currently observed blobs are merely 
among the brightest members of their class at z = 3. Our derived luminosity function goes roughly 
as dM(> L)/dL oc L"^'^ at S99's luminosity threshold. If the sky background dominates the pixel 
noise, detecting much fainter objects at 2; = 3 will be difficult. However, our luminosity evolution 
suggests the surface brightnesses of the blobs may peak at 2; ~ 1. Even at z = 2, we expect ~ 6 
times as many blobs at S99's limiting surface brightness as at 2; = 3 (these numbers are sensitive 
to cosmology). Ground-based observations at the lowest achievable redshift, z 2, may be the 
optimal method for finding these blobs until the launch of NGST. 

The results we find have applications beyond the Lya blobs. For example, semi-analytic models 
of galaxy formation generally assume that gas in galaxies should be cooling principally at the virial 
temperature of the galaxy, and hence be emitted in X-rays. In the picture of White & Frenk (1991), 
the gas in a halo shock heats to the virial temperature, then condenses onto the central galaxy from 
the inside out as it cools. The rate of cooling is given by the minimum of the dynamical gas infall 
rate and the growth of the mass within the cooling radius, defined to be the radius where the cooling 
time equals the cosmic time. The emitted cooling is taken to have a near isothermal spectrum at 
the virial temperature. 
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A problem with this picture is that the expected X-ray emission from individual galaxies has 
not been detected with ROSAT (Benson et al. 1999). One reason for this failure may be that the 
cooling radiation is not emitted primarily at the virial temperature but instead covers a large range 
of temperatures, with 10^ K predominant. This in turn suggests that the assembly of gas to form 
galaxies is a more gentle process than in the White & Prenk (1991) picture. Once gas has collected 
into cool dense objects, it seems that it is difficult to force it out of that state, and dissipation 
mostly takes place through efficient atomic lines. Dominance of low temperature cooling does not 
necessarily require that galaxy assembly proceed by mergers of discrete, cold objects. If the gas 
encounters a series of weak shocks as it falls into a galaxy instead of a single strong shock, it may 
be able to radiate its internal energy as quickly as it acquires it, never reaching high temperature. 

During preparation of our paper, a preprint by Haiman et al. (2000) appeared proposing a very 
similar model for the S99 blobs. We note some of the similarities and differences of the two models 
here. Haiman et al. (2000) base their analysis on a semi-analytic calculation, using a merger tree 
formalism and a spherical collapse model. In contrast to previous semi-analytic calculations that 
assume the energy emerges at the virial temperature, Haiman et al. (2000) argue that cooling will 
be so efficient that most energy will be radiated at ~ 10^ K. In general, their picture is in reasonable 
agreement with ours, but it differs in quantitative details. We find lower amounts of Lya radiation 
by a factor 2-4, partly because they ignore H I 2-photon emission and bremsstrahlung as coolants. 
Their treatment of radiative transfer seems to overemphasize the frequency shift due to resonant 
scattering and underestimate the spatial scattering. Because they assume the Lya radiation is 
produced out to the virial radius, they do not consider this a problem. Our emission is somewhat 
more centrally concentrated. They assume monolithic collapse, so that for a given galaxy there is 
first a cooling stage and then a star formation stage. In contrast, we find that cooling and star 
formation happen simultaneously in massive galaxies. 

Finally, Haiman et al. (2000) find that the efficiency of star formation in producing Lya 
radiation is only twice that due to gravitational cooling. In our simulation this ratio is about 20. 
The difference comes in part from their channeling all the cooling radiation through Lya and in 
part from their assumption about the importance of feedback. They assume that only 10% of the 
gas forms stars. Since at the present day most galaxies are dominated by stars rather than gas, 
this implies the remaining gas must be blown out again by galactic winds. In our simulations, the 
efficiency of forming stars out of eligible gas is only 10% at each timestep, but the gas is not blown 
out, and eventually most of the cooled gas is turned into stars. The effectiveness of feedback in 
removing gas from galaxies is currently a major question in astrophysics, and the comparison of 
cooling radiation and stellar emission may eventually help to constrain the answer. 

Further theoretical work is needed to test the robustness of our predictions over a wider 
dynamic range of numerical parameters. Further observational work is needed to test whether the 
predicted Lya emission from young galaxies exists in the real universe. This emission, perhaps 
already observed in the form of the Lya blobs, represents a novel form of radiation from galaxies 
that potentially offers a direct view of the process of galaxy formation. 
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Fig. 1. — Cumulative distribution of the cooling radiation as a function of the temperature of 
the emitting gas, for the L144 simulation at z = 3 (solid) and z = (dashed). The supernova 
contributions to the cooling are omitted. The dotted line shows the distribution expected for gas 
that cools from an initial temperature of 10^ K. Since the gas is fully ionized in H and He down to 
^-^8 X 10*^ K, this curve is approximately T/(10^ K) over most of its range of significance. 
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Fig. 2. — Luminosity functions based on the different sources of energy input to the gas. The 
gravitational cooling luminosity function of groups in the L144 simulation at z = 3 is plotted as 
the dot-dashed line. The contributions from supernova heating and stellar Lyman continuum are 
plotted as the dotted and dashed lines respectively. The sum of gravitational and supernova inputs 
are plotted as the solid curve. The turnover in the curves at a density of 3 x 10~^ Mpc^^ is an 
artifact of limited resolution. All luminosity functions in this paper are plotted using comoving 
densities. 
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Fig. 3. — Evolution of the cumulative Lya cooling luminosity function at several redshifts, for the 
L144 simulation. Only the contribution from gravitational sources is included. 
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Fig. 4. — An example of the Lya radiation resulting from gravitational cooling (left panel) and 
stellar photoionization (right panel) in one of our groups at z = 3. The example here is taken 
from our highest resolution simulation LI 1/64. The emission is shown by the gray scale images; 
the intensity scale is marked in terms of log]^Q[/ij^Q,/(ergs~"'^ cm~^ asec~^)]. The neutral hydrogen 
column density is also shown at contours of logj^Q(A^Hi/ cm~^) = 18, 19, 20, 21, and 22. 
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(gas) (proper, h ^ kpc) 

Fig. 5. — Radial rms extent of the gas in the cooHng groups at z = 3, weighted by several different 
factors. The x-axis is weighted by the gas mass. On the y-axis, the triangles are weighted by 
the Lya emission (considering only the "gravitational" cooling). The squares are weighted by the 
star formation rate, which is proportional to the stellar ionizing emission. The dotted line shows 
the gravitational softening length. The points would lie on the solid line if the emissivity were 
proportional to the gas density. Only groups with Lya cooling of more than 10^^ ergs~^ are shown; 
groups with Lya cooling of more than 5 x 10^^ ergs~^ are shown surrounded by outlines. 
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Fig. 6. — Cooling as a function of the "galactic mass", at z = 3. This is defined as the mass of the 
gas with temperature T < 3 x 10^ K and overdensity pg/pg > 1000, plus the mass in stars. The 
cooling includes only the gravitational Lya cooling. The galaxies in the L144 simulation and the 
Lll/64 simulation are plotted as squares and triangles respectively. 
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Fig. 7. — The total gravitational cooling radiation as a function of the star formation rate, at 
z = 3. The galaxies in the L144 simulation and the Lll/64 simulation are plotted as squares and 
triangles respectively. The Lyman continuum luminosity resulting from the star formation is shown 
along the top axis. The supernova heating resulting from star formation is shown by the dashed 
line; galaxies lying on this line have equal gravitational and supernova contributions to the cooling 
radiation, and gravitational cooling dominates for points above the line. Galaxies on the dotted 
line would have gravitational cooling radiation equal to the ionizing radiation output of their young 
stars. Galaxies with no star formation are plotted at the extreme left. 
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Fig. 8. — Total gravitational cooling vs. potential energy, at z = 3. The galaxies in the L144 sim- 
ulation and the Lll/64 simulation are plotted as squares and triangles respectively. The potential 

energy here is simply calculated from U = Mjl/^^^, where Vdrc is the measured circular velocity 
and Mb is the baryonic mass. The dotted line shows the potential energy divided by the cosmic 
time. 
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Fig. 9. Cooling luminosity in Lya at two resolutions at z = 3. Galaxies from the lower resolution 
Lll/32 simulation, which has the same resolution as L144, are plotted versus the same galaxies 
from the higher resolution, Lll/64 simulation. 



